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ABSTRACT 


Analytical and numerical models are developed to predict 
the longshore sand transport distributions across the surf 
zone. The models, which also predict the root mean square 
wave height, Has! and the longshore current, V, are compared 
with field data acquired from Leadbetter Beach, Santa Barbara, 
California, during the intensive storm period in February, 1980. 
The breaker coefficient, B, and the bed shear stress coefficient, 
Cer when equal to 1.1 and 0.005, were found to give the best 
agreement between the predicted aan and V with the field 
measurements. The model for a plane sloping beach predicts 
the maximum sand transport at offshore distance X equal to 
0.8X,, where Xp 1s the mean breakerline location, which agrees 
with Komar (1977a). The Longshore sand transport formula, 
suggested by CERC (1977), is used to calibrate the model which 
requires the empirical transport coefficient Bre EOE. Ole 

The sand transport cross-shore distribution predicted by 
the model does not agree well with the corresponding trans- 
port inferred from the field measurements. Reasons offered 
to explain the differences are (1) that the model does not 
include the swash zone where the maximum transport is found 
and (2) that beach profile information beyond a distance of 
100 m offshore was incomplete. The results indicate the 
importance of including the swash zone to describe the effect 


of the long waves in predicting the sediment transport. 
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io INGRODUCTION 


Water waves are one of the principal causes of shoreline 
changes. When waves break along the shore, they release their 
energy and momentum and give rise to a longshore current. 

The longshore current, along with the stirring action of the 
waves, is the primary mechanism for longshore sediment trans- 
port. The longshore sand transport rate is an essential fac- 
tor determining erosion or accretion along a coast. If the 

longshore sand transport rate can be accurately estimated, a 
quantitative picture of shoreline evolution can be evaluated, 
including changes of the shoreline due to marine structures. 

In the present study, analytical and numerical models are 
developed based on a longshore current model for random waves 
by Thornton and Guza (1983) and a sediment transport formula- 
tion by Thornton (1973) to predict the cross-shore sediment 
transport distribution and to compute the total volume of sand 
transport rate. The model is compared with the field data 
acquired from Leadbetter Beach, Santa Barbara, California 
(January-February, 1980). 

The experiment at Santa Barbara was conducted as part of 
the Nearshore Sediment Transport Study (NSTS). The objectives 
of NSTS are to develop improved techniques and engineering 
formulae to predict sediment transport on beaches having 
Straight and parallel contours by utilizing field measurements 


of waves and currents at several nearshore Sites. 


alia 





During the Santa Barbara experiment, beach profiles were 
measured daily along with extensive measurements of waves and 
currents. A severe storm attacked the site of the experiment 
during 16-20 Feb., 1980 causing significant changes in the 
shoreline configuration and beach profiles. Despite the 
destruction of most of the instruments during the storm, 

a very substantial and valuable record of waves and currents 
was acquired which makes the Santa Barbara data very unique and 
Significant for nearshore sand transport studies. 

The present work emphasizes the data acquired on 17 and 
18 February, 1980. A strong longshore current was observed 
due to storm waveSapproaching the beach at moderately large 
angles (6° at the breaker line) causing significant erosion 
in the beach profile. The daily beach profiles during the 
experiment (January-February) show no well defined bars or 
troughs and essentially have a stable point 95 m offshore, 
which did not appear affected by erosion or accretion. Thus, 
it is hypothesized that the changes in the beach profile are 
due to sand transport in the longshore direction alone, and 
the on-offshore transport is negligible. Furthermore, the 
longshore sand transport distribution across the suri zone, 
predicted by the model, is assumed to have a distribution 
Similar to the net longshore sand transport calculated from 


the profile changes. 
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II. LONGSHORE SEDIMENT TRANSPORT 


A. GENERAL 

When waves approach the shoreline at an oblique angle and 
break, they generate a longshore current parallel to the 
shoreline. This current in turn interacts with the waves to 
produce a longshore sand transport. The movement of sand 
particles is usually divided into two, not always distinct, 
modes of transport: bed load and suspended load transport. 
The bed load is the result of sand grains being rolled along 
the bottom by the shear of water moving above the sediment bed 
and 1S maintained by grain-to-grain contact. Bed load is 
usually considered to be confined to within 10 cm of the bottom. 
In the suspended load mode, the generally finer grains are 
transported by currents after they have been lifted from the 
bed by turbulent action of the wave-induced vertical motion. 
The total sand transport is the sum of the contribution by 
bed and suspended load transport. The rate at which the total 
volume of sand is transported parallel to the shoreline is 
termed "longshore sand transport: Q,". The two most success- 
ful and widely used approaches for estimating the longshore 
sand transport are an empirical wave power model and 


Bagnold's energetics model. 


B. WAVE POWER MODEL 
In the wave power model, the rate of total longshore 


transport Q, is simply assumed proportional to the longshore 
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component of wave power (P,) at the breaker line to some power 


n (Watts, 1953a). 
Qo =e ee (221) 


where Ky and n are coefficients to be determined empirically. 
Assuming straight and parallel depth contours, the onshore energy 
flux in the direction of wave propagation per unit length of 

the beach is EC, COs a, where E is the energy density of the 
waves, Os is the group velocity of the incoming waves, and 

a is the angle of the wave crests relative to the shoreline. 
Then, the longshore component of wave energy flux per unit 


length of the beach is given by: 


Pp = VEC a) js SEK Che, SEE Oh a ) 
where the subscript b refers to the breaker line. Applying 
Hinear theory, 
P = eae COS aS 10 cl (2635) 
Q 8 DEG b b : 
where op is the water density, g is the gravitational con- 
stant and H is the wave height. 
Watts (1953a) obtained the first field measurements by 
which the sand transport rate could be related to the local 


wave characteristics. Caldwell (1956) combined additional 


ie 





data with those of Watts and found the coefficient (n) in 
weet) co be less) than unity, and with a different proportion- 
ality coefficient (Ky) than Watts. 

Inman and Bagnold (1963) interpreted field and laboratory 
data and obtained a linear relationship between Q, and Poy 
i1.e€., n equal to unity, which takes into consideration the 


immersed weight of beach material. 


K 

= P (2.4 
°c“ TI=ay(ep-e) Pe | 
where - is the density of sand transport, ais the porosity 
of beach sediment (taken as 0.40) and K is a dimensionless 
Proportionality coefficient. 

Komar and Inman (1970) utilized sand tracers to measure 
longshore sand transport rates. Using this new data combined 
with data from Watts (1953a), Caldwell (1956) and others, they 
obtained 0.77 for the proportionality coefficient K by relating 
the immersed weight transport rate I, to Py, 

I = 0.77 P (2725) 
3 3 
For quartz-density sand, De bs 2eoo ieee kGg/m 5 K = 0.77, 


equation (2.4) becomes, 


Q, ro. OP 276) 


16 





es 


watts/meter. Komar and Inman use the root mean square wave 


where Qy in (2.6) is measured in units of m>/day and P 


height in (2.5) and (2.6) instead of the significant wave 
height, suggested by CERC (1977), to evaluate the Py term. 
The relationships described in equations (2.5) and (2.6) 
are almost purely empirical with no real consideration of 
the sand transport mechanism. In addition, the coefficient 
0.77 in equation (2.5) has a high degree of uncertainty 
Longuet-Higgins (1972), in attempting to rationalize (2.2), 
points out that since the energy flux VES is a vector rather 
than a second order tensor, the longshore component of wave 


power (P,) can be written equivalently as the product of two 


physically meaningful quantities: 


ee: 
= sh 4 — 
Po [E a cos a sin a]Cc Sy C C2) 


where Soy is the longshore component of the radiation stress 
and C is the phase speed of the waves. The radiation stress 


(S,. ) is shown by Longuet-Higgins (1972), among others, to be 


iy 
the principal driving force for longshore currents, and thus 


would be expected to be important in describing longshore 
sediment transport. Vitale (1981) compared the radiation 
stress ee 


) and the wave power (P,) as predictors of the 


Y R 
longshore sediment transport rate by conducting three-dimensional 
movable-bed laboratory tests and he found no major difference 

in choosing S sey over Py to predict the longshore sediment 


meansport. 
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C. BAGNOLD'S ENERGETICS MODEL 

Bagnold (1963) relates the rate of immersed weight trans- 
wort (I, ) to the work done by waves and current. The mechanics 
of sand particle movement are described as a back-and-forth 
motion under the action of the wave orbital motion with 
essentially no net transport. Wave energy is expended in 
Supporting and suspending the moving sand above the bed. 
Once the sediments are suspended, the presence of a unidirec- 
tional current superimposed on the to-and-fro motion can then 
produce a net drift of sediment. Bagnold (1963) derived the 


relationship: 


"9 
1, = K Ww > (2528)) 
O 


where i, 1s sediment transport per unit width in the direc~ 

tion @ determined by the unidirectional current Ug, Ww alisy elev 
available wave power supporting the sediments above the bot- 
tom, U, is the orbital velocity of wave motion, and K' is a 

dimensionless coefficient. 


Inman and Bagnold (1963) specified the various parameters 


of waves at the breaker line and obtained: 
(25. 95 


where u, is the maximum horizontal orbital velocity of the 


waves. Equation (2.9) can be viewed as a general relationship 
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because it does not specify the cause of the longshore cur- 
rent Vpi it can be either tidal, currents of a cell circula- 
tion, wind generated, or due to oblique wave approach. Komar 
and Inman (1970) utilized their littoral drift measurements 
to test equation (2.9) and required K' to be equal to 0.28 

to agree with their data. In terms of Qy and for quartz 


density sand, equation (2.9) is given by 


Yb 
Qe =e} SENS oie Tad 22 erle)) 


= 


where Q, 1s in m>/day and the wave and current parameters 
are in mKs units (watts/meter) (Komar, 1983). Wang and Chang 
(1979) found good agreement between their measurements along 
the bayshore of a barrier island in the Gulf of Mexico, and 
equation (2.9) for a coefficient K' = 0.18, where the long- 
shore currents generally are due to local winds and tides. 
Bowen (1981) applied Bagnold's equation to the problem of 
on-offshore sediment transport on beaches and then studied 
the special case of a normal incident wave and steady on- 
offshore currents (no longshore currents). The waves are 
described using Stokes second order wave theory and the mean 
on-offshore current using Longuet-Higgin's (1953) bottom 
streaming solutions. The model predicts an equilibrium 
beach profile as a function of the incident wave character- 
istics. Bowen's results support observations that steep 
beaches are generally coarse grained, and short period waves 


Cause a beach to erode. 


aR, 








Pagure (2.1). 


Schematic of the Bagnold model for sand 
transport wherein the orbital velocity 
UQmeoe co Heewaves places the Sand in 
NCct Om menmcncmeurrene Us provides a 
net transport of sand, I). 
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Bailard and Inman (1981) independently derived a very 
Similar model to Bowen's (1981) to describe the instantaneous 
transport in the nearshore zone, including both on-offshore 
movements as well as longshore movement. An important part 
of the study is the effect of local bottom slope on sand 
transport rate, as this permits the analysis of the formation 
of the longshore bars and other inshore topography. Bailard 
(1982), extended the model to predict the total load trans- 
port for time-varying flow over an arbitrarily sloping planar 
bed. The model predicts the local equilibrium beach slope 
as well as the local near bottom sediment transport rate as 
a function of the near bottom velocity. Surf zone conditions 
are described as steady longshore and on-offshore currents 
in addition to a wave induced current having a local wave 
angle. 

A number of other models have been proposed. For instance, 
Madsen and Grant (1976) adapt the Einstein-Brown sediment 
transport equations to time-varying sand movements under com- 
bined waves and unidirectional currents. The bottom stress 
(tT) under the combined oscillatory motions of the waves and 
any superimposed currents, and the resulting sand transport, 
vary with time. The model is relatively complicated and re- 
quires parameters for which no data are available. The 
Success of the wave power formulation and Bagnold's energetics 
model is due to their relative simplicity requiring only one 


or two parameters specified from the data. 


Ze 





ime general, the predictions of the littoral drift have a 
mgm deqree Of uncertainty. This partly results from our 
basic inability to make accurate measurements of sand trans- 
port on beaches. However, the two approaches of equations 
(2.5) and (2.9) can be used to give at least a qualitative 


Mescription of littoral drift. 


tee DLESTRIBUTION OF LONGSHORE TRANSPORT ACROSS THE SURF ZONE 
Bagnold's (1963) approach has been also utilized to pre- 

dict the distribution of longshore sand transport. Thornton 

(1973) analyzed sand transport inside and outside the surf 

zone where the volume transport rate per unit width (q.) 1s 

given by: 

Be a WES 


g(l - = On 
S 


e272 nl 


1Q | 
n 
Il 


aX 


where ae 1S an empirical dimensionless coefficient. Thornton 
(1973) obtained field data at Fernandina Beach, Florida, to 
evaluate B_. in equation (2.11). The transport rate (qo) was 
measured using a series of bed load traps placed in a line 
across the nearshore. The results indicate that the maximum 
longshore transport is in the breaker zone with Be = 0.08 
and stronger transport occurs over the bars rather than in 
mee troughs (Fig. 2.2). 

Komar (1976b,1977a) similarly utilized Bagnold's (1963) 


meael £O analyze the distribution of the littoral drift in 


2 
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The distribution of the longshore sand 
transport measured by Thornton (1973) 
at Fernandina Beach, Florida, by bed 
load traps in the positions as shown. 
The dashed lines are based on equation 


Zeer WithB-. =.0.08. [from Thornton 
(1973) ] = 
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terms of (I, ). The local immersed-weight of sand transport 


rate per unit of width inside the surf zone (1,) is given 





Dy : 
ie 2 
i(x) = 7 (On Ce) OY DU SG AU UR) (242) 
where C,. is the drag coefficient of oscillatory wave motions 


is 
and h and v are the local water depth and longshore current, respec- 
tively; both are functions of the distance xoffshore from the shore- 
line. Ky is a dimensionless proportionality factor deter- 


mined by calibrating the model equation (2.12) to yield I, 
calculated from equation (2.5) where the total immersed- 
weight transport can be evaluated by integrating i(x) across 


the surf zone: 


where Xp Smee uSUint 2ZOMe wldtn aS SnOwn In Fig. (2.3). 
Komar (1977a) calculated the distribution of both the long- 
shore current and the longshore sand transport, Fig. (2.4), 
based on the longshore current analysis by Longuet-Higgins 
(1970). The distribution shows the maximum longshore current 
at approximately (X/ xX, = 0.65) and the maximum longshore sand 
transport at approximately (X/X) elo 0))) 2 

Sawaragi and Deguchi (1979) utilized circular bed load 


traps which enable them to measure the on-offshore sand 
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meguee  (2.3)5 Schematic of defining the location and 
the water depth at the breaker line 
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Figure (2.4). Longshore current and sand transport 
distributions determined by Komar 
e774). 1 fsom, Komar (19774) |] 
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movements as well as the longshore transport distribution 
across the surf zone. The results showed sediment transport 
to be relatively independent of wave steepness Ho/h: but 
dependent on the sand grain size. The maximum sand trans- 
Meee is shown (Fig. 2.5) to occur at approximately 
(X/X, = 0.6) and maximum longshore current at (X/Xp = 0.4). 

Bailard (1982) extended his earlier work in which he used 
second order Stokes theory to describe the wave motion. He 
included the Ostenddorf and Madsen (1979) longshore current 
model to predict the distribution of sediment transport rate 
for different ratios of the velocity of the flow to the sedi- 
Memienrall velocity (Fig. 2.6). He finds the maximum long- 
shore current at approximately (x/X, = 0.6) with the maximum 
longshore sand transport at (X/X) = Oige es 

Results from different studies examining the longshore 
current and longshore sand transport distributions are given 


mae Table I. 
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Figure (2.5). Laboratory measurements of the distribu- 
tion of the longshore current and the sand 
transport rate showing the positions of 
the maximum sand transport and maximum 
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tito DESCRIPTION OF LEADBETTER BEACH EXPERIMENT, 
SANTA BARBARA, CALIFORNIA 
A. GENERAL 

Leadbetter Beach, Santa Barbara, California was chosen as 
the second NSTS experimental site, after Torrey Pines Beach, 
California, primarily because it is adjacent to a total sand 
trap forming a spit off the Santa Barbara breakwater, and 
because it has relatively straight and parallel nearshore depth 
contours. The NSTS Santa Barbara experiment began in October, 
1979 with a 14-month sediment trap study measuring sand 
accumulation, coupled with a one-month intensive experiment 
during 27 January through 25 February 1980. 

Exceptional storms were encountered towards the latter 
part of the intensive experiment. During the period 16-20 
February 1980, the weather conditions deteriorated greatly 
approaching the 50-year storm. Significant wave heights in 
excess of 2 meters and longshore current speeds of 150 cm/sec 
were recorded. These conditions coupled with spring tides, 
high winds and wave set-up, eroded the foreshore and back 
beach areas approximately 2 meters in elevation. These 
phenomena coupled with high quality data, make the Santa 
Barbara experiment very unique and significant for near-shore 
transport studies (Gable, 1981). This study emphasizes the 


very intensive storm period of 17-18 February, 1980. 
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B. COASTAL GEOMORPHOLOGY 

The shoreline between Point Conception and Santa Barbara 
has the unusual (for the California Coast) east-west orienta- 
tion on a predominantly north-south coast. The coastal areas 
of Santa Barbara are composed of marine sedimentary rock and 
Monterey Shale-bluffs 30 m high fronted by sand and cobble 
beaches varying from 0 to 30 m in width (Corps of Engineers, 
1970). The major sources of sand for the Santa Barbara 
beaches are the Santa Maria and Santa Ynez rivers (see Fig. 
3.1). Some sand comes from the erosion of the coastal cliffs 
(Trask, 1952). Bascom (1951) suggests some contribution by 
wind blown sand, but neither of these two sources are rela- 
tively important. Textural analysis of the sands indicates 
the beaches are composed of fine to medium grain size and are 


relatively well sorted with a mean diameter of 0.236 mm. 


C. WEATHER 

The climate in the Santa Barbara area is classified as a 
Mediterranean type and characterized by warm, dry summers and 
mild, wet winters. The summers are mostly dry because a 
semi-permanent high pressure area covers the Eastern North 
Pacific Ocean, which deflects eastward storms to the north. 
During winters, this Pacific high migrates southward and 
weakens, allowing occasional frontal systems to move through 
Southern California. 

The prevailing winds for the West Coast are from the 


northwest due to the semi-permanent Pacific high pressure and 
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blow generally parallel to the coastline. However, due to 
the orientation of the coastline at Santa Barbara and the 
Pmeoceccion from Point Conception, the winds are generally 
Peeiemanidevaridele rrom the Southerly direction. The aver- 
age wind speed is less than 3 m/sec and the typical pattern 
meethat Of a breeze. 

Most extratropical winter storms that affect the Santa 
Barbara area move southeastward from the Northeast Pacific; 
on the average, 1 or 2 such storms influence the Santa 
Barbara channel region every year. 

In February 1980, the Eastern Pacific high pressure sys- 
tem was weaker than normal and displaced southeast of its 
normal position. At the same time, the low pressure center 
was also displaced south of its normal position. The surface 
pressures were as much as 18 millibars below normal east of 
the low center. The Eastern Pacific high pressure system 
was almost non-existent, which shifted the storm track suffi- 
ciently to bring storms onshore in Central and Southern Cali- 
memmta (Miller, 1980). 

Two series of frontal systems hit the study site. A 
memong SOucheasterly occurred from 13 to 16 pepeaeey and a 
very intensive southwesterly from 17 through 21 February. 
Because these storms were southerly, they produced large 
westerly ocean swells and waves which moved into the Santa 
Barbara channel. The relatively warm water was a major factor 


which intensified the storm and brought much heavier rain 


2) 





than normal. Precipitation was 13.2 cm above normal at this 
time of the year. All these factors, coupled with extreme 


high tides and wave height, caused extreme beach erosion. 


D. WAVES AND WIND CLIMATE 

The wave climate at the experimental site is either locally 
generated inside the Channel Islands or propagated from the 
open ocean. Normally, the predominant waves are swell from 
the west-southwest and the west (240°-270° true) entering 
through a narrow window between San Miguel Island and Point 
Conception with an average wave height of 91 cm and an average 
wave period of 12 seconds (Weigel, 1950). The waves must 
refract almost 90° to reach Leadbetter Beach (Fig. 3.2). As 
a consequence, the swell waves reaching Santa Barbara are 
greatly reduced in height and arrive at large angles due to 
the severe refraction (O'Brien, 1950). 

As a result of the origin of the waves, the offshore 
islands, and the refraction of the waves, the wave pattern 
in Santa Barbara is nearly the same for all times of the year, 
with the exception of the waves from the occasional southeast 
storms. The swell from the open ocean, refracting and break- 
ing at high angle of incidence, create a strong easterly 
longshore current. The occasional southeasterly storms also 
create reversals where the longshore current is to the 
west. 

The local wind waves from the south are usually insignifi- 


cant with only choppy seas and small waves. The local storms 
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from the southeast (125°-145° true) have a maximum fetch of 
145 Km and can propagate waves toward the Santa Barbara Coast 
with significant wave height ranging up to 488 cm (Weigel, 
1950). During the storm period (16-20 February, 1980), waves 
generated thousands of kilometers away with unlimited fetch 
were funneled through the wave window of (240°-270° true). 
The wave period was 14-16 sec with significant wave height 
approaching 200 cm from the southwest and west. The signifi- 
cant wave height was 190 cm from the southeast on February 16 


and 185 cm from the west on February 20, 1980. 


E. EXPERIMENT 

The start of the experiment was timed such that a strong 
negative low tide occurred during the daylight hours allowing 
nearshore instruments to be installed on a dry beach. The 
tides are semidiurnal with a mean range of 110 cm referred 
to mean lower low water. The tides ranged from -10 cm to 
+200 cm (MLLW) during the intensive experiment. 

A plan view of Leadbetter Beach (Fig. 3.3) shows the loca- 
tion of the instruments, the experiment baseline and the five 
range lines. It is noted that the five range lines are not 
equally spaced. 

Series of current meters, pressure sensors and run-up 
meters were used to measure the wave characteristics during 
shoaling, breaking, and run-up. The location of the two wave 


slope arrays is shown in (Fig. 3.4). 
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l. Incident Wave Field Measurements 

Two, 6 meter square arrays, composed of four Kulite 
solid state pressure sensors were used to measure an un- 
biased estimate of the spectrum of the radiation stress, S.. 
The square arrays are less efficient than a linear array in 
estimating a complete directional spectrum but the narrow 
wave window provided by the offshore islands of Santa Barbara 
Beach produces an acceptable wave climate estimate. The 
location and orientation of the two slope arrays are shown 
in Fig. (3.4). The water depth at both arrays 1s approximately 
9 m and the sensors were mounted 120 cm above the bottom. 

2. Profile Measurements 

Profile measurements at Santa Barbara consisted of 
daily wading profiles of the beach face and shallow nearshore, 
coupled with offshore surveys to a depth of 4 m using a trac- 
tor and boat survey for pre- and post-experiment measurements 
to 12 m depth. Five range lines, extending essentially per- 
pendicular to the beach and covering 200 m alongshore, were 
established for the daily wading profile (Fig. 3.3). The off- 
shore distances along range lines are referenced to the distance 
from the experimental baseline (y-axis). 

The wading profiles were done using conventional rod 
and level survey techniques on a daily basis. Range lines 
were defined by two survey points marked by flags, which a 
rodman aligned before each survey reading. Rod stations were 
measured using a specially constructed plastic-coated steel 


survey line marked at 5 m intervals. An engineer's level and 
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rod provided vertical control. The profile survey was ter- 
Minated when the water became too deep for the rodman, or 
the breaking waves made it impossible to plumb the rod. These 
Surveys were completed at the lowest possible tides occurring 
during daylight hours. Each survey was closed to the original 
benchmark on land to verify survey accuracy. Data were 
entered on site by hand into a Tektronix computer terminal 
within a short time after each survey. Profiles were plotted 
immediately to visually search for errors which were due to 
either data entry or a bad survey point which could be cor- 
rected by redoing the field survey immediately. 

Generally, small scale changes occurred in the profile 
during most of the intensive month-long experiment (Fig. 3.5), 
except during the storm where the beach was cut vertically 
2.5 m and horizontally about 60 m withinalmosta four day period. 
Although it was impossible to maintain the nearshore instru- 
mentation during the storm, the deep water instruments con- 
tinued to work, providing an extremely valuable data set on 


beach erosion. 
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IV. FIELD DATA ANALYSIS 


A. BEACH PROFILES AND CONSERVATION OF SAND TRANSPORT EQUATION 
Five beach normal ranges at Leadbetter Beach were surveyed 

on a daily basis starting from the established baseline on 

the beach to a distance about 100 m offshore. Standard rod 

and level technique described in (III.E.b) were used. Data 

from the boat-fathometer survey conducted on 25 February, 1980 


for the five ranges were used to extend each profile out to 


1 
59) - 


Linear interpolation was applied twice, first to calculate 


the shallow water limit at x = 180 m where é < 


the missing data points, and second to define five equally 
Spaced ranges (50 m apart). The five ranges were averaged 
for both 17th and 18th of February to produce a daily mean 
profile. Then, a “mean beach" profile was defined by averaging 
these two daily profiles. Fig. (4.1) shows the mean beach 
profile, relative to mean sea level (MSL) and indicates no 
well developed bar structures or trough along the profile. 
Based on the difference in elevations between the 17th 
and 18th February surveys at each location, the net longshore 
transport distribution is calculated using the conservation 
analysis of sand transport described later in this section. 
A "mean net longshore transport", q_.,, distribution was ob- 


tained by averaging the net transport distributions for the 


five ranges. The "mean net sand transport" distribution and 
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the standard deviation calculated at each point are shown in 
Pig. 4.2 indicating a small variance with acceptable scatter 


around the gq lal Geral abla sveygle 
mean 


The Gee distribution goes to zero after a distance of 


100 m offshore as shown in Fig. 4.2 because the wading pro- 
files on the 17th and 18th are extended after this distance 
by the same fathometer data taken on 25 February. 

Profile changes between sequential measurements can be 
approximated by a parallelopiped with sides Ax, Ay and the time 
varying depth (Ah). Considering the equilibrium of the 


parallelopiped element (Fig. 4.3), the net sand trnasport rate 
og 
x 


Sey = 5x f xy and the net sand 


in the x-direction is q 


Tet 
transport rate in the y-direction is g Ax = yhyax. 


Vane te 


The total net sand tranSport rate (g hes as 


Pie ies dy nee 


balanced by the rate of the net volume changes over the ele- 


ment giving 


= Ah | 
Tenet SY + Gy net 4* = (l-a) 7p Axdy (4.1) 
where a is the sand porosity and is taken as 0.40. Equation 


mel) Can be written as 


dq dq 
See ay on 
oe + a (l-a) = (422) 


| Q> 


Since a strong longshore current was observed during 17 


em@cm'G February, (see Fig. (6.3)), a significant longshore 
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Prgure (4.3). Conservation of sand transport with 
changes in the beach profile. 
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sand transport was expected. Therefore, it is hypothesized 
that the longshore sand transport is much larger than the on- 
offshore sand transport, i.e., ah >> dy. The profile changes 
in Fig. 3.6 do not show any well-defined bar or trough 


structures which suggests the on-offshore sand transport is 


3q 
negligible (i.e., == = 0). Thus, equation (4.2) becomes 
oq 
nanan ee 
= (l-a) TE C4733) 


Equation (4.3) states simply that the changes in the beach 
profile with respect to time, areproportional to the net 
longshore sand transport. The transport rate/unit width of 
on-offshore distance, ay can be calculated since water depth 
(h) is measured as a function of time (t). The boundary con- 
ditions are needed to solve (4.3). For the offshore boundary 


condition it is assumed 


where Ae is the water depth at the shallow water limit, 


h and L is the wavelength. Given the dominant wave 


je 
O 20’ 
period of 14 sec (f — On O 7 Os) Hz). hy was Calculated to be 
5.0 mat a distance approximately 180 m from the baseline. 
The nearshore boundary condition requires the sand transport 
to go to zero at the shoreline. It appears from Fig. 4.2 


that this condition is satisfied from the data at the baseline 


at x = 0. The net sediment transport rate is calculated by 
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integrating the net cross-shore distribution which is assumed 
to be similar to the longshore sand transport distribution. 
For discussion purposes, two regions within the surf zone 
are defined. A swash region is defined from the baseline to 
the intersection of HWL with the profile, i.e., (from x = 0 
to x - 30 m); a seaward region is defined from x = 30 m 
to x - 180 m. The swash region indicates significant 
erosion on the beach face (Fig. 4.2), which will not be considered 
in this study because the model does not include sediment 
transport on the beach face. The seaward region shows the 
eroded profile with a small accreted area around x = 90 mM, 
which appears with a negative Sign. 
Bee RADIATION STRESS, WAVE HEIGHT AND DIRECTION AT OUTER 
BOUNDARY 
Measurements from a 610 cm square, four pressure sensor 
array in 8-m water depth are used to produce an unbiased 
estimate of the spectrum of radiation stress. The pressure 
arrays worked throughout the experiment including the storm 
period. The array measurements are also used to calculate 
the height and the directional characteristics of the waves. 
The radiation stress (Sy) 1S given by (Longuet-Higgins, 


Hs 6 2 ) 


oe - pu (z,t)v(z,t) dz CATA 


where u, and v., are the crossshore and alongshore components 
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of wave velocity. The velocities are calculated from the 
pressure gradients measured at the array applying linear wave 


theory spectral transformations. For a particular component 


= = ri cosh k(h + z), dp nes 
u,(£,2,¢) = ATEOGIN aN a z ] ax, = = See (4.5) 


where the Za is the measured depth of the pressure sensors, 
W is the radial frequency and j refers to the two horizontal 
components. The term between the brackets is the spectral 


transfer function H(f) relating the measured pressure gradients 





to inferred velocity components. Substituting equation (4.5) 
into (4.4), using the meaSured pressure head = and integrating 
Sen(f) = ee) G Le) (4.6) 
ey 2u'k(cosh*k(h+z)) 3B, of 


where Cas apf) is the co-spectrum between the measured 


ax’ Dy 
pressure gradient components. 


Wave heights were determined from the surface elevation 
records using the "Zero-up-crossing method" in which the wave 
height is defined as the difference between the maximum and 
Minimum occurring between two consecutive zZero-up-crossings. 
surface elevation records were first linearly detrended to 
exclude the effects of the rising and falling of the tides 
and then high-pass filtered with a cut-off frequency of 0.05 
Hz to exclude surf beat, followed by a low-pass filter with 


a high frequency cut-off (0.3 Hz). 
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The filtering was accomplished by Fourier transforming the 
signals, zeroing the Fourier amplitude coefficients in 
the filtered frequencies, and inverse transforming the 
complex spectrum to obtain the filtered time series. The 
entire 68-minute-record was transformed at one time to mini- 
mize the end effects which result in spectral leakage, and 
to obtain maximum resolution giving very sharp roll-off at 
the filter cut-offs. 

Due to the large number of current meters, current data 
are utilized applying linear theory to infer wave heights 
from the velocities (Guza and Thornton, 1980). Surface ele- 
vation time series were obtained from current measurements 
using the complex Fourier spectra of the horizontal velocity 
component U(f), V(f). The complex surface elevation spectrum, 
X(£), is calculated applying the linear wave transfer function, 


H (£) 
X(£) = H(£) + V(£) (4.7) 


Then, the complex surface elevation spectrum is inverse 
transformed to obtain the surface elevation time series from 
which the wave height distribution is calculated. Surface 
elevations were also inferred from pressure signals in a simi- 
lar manner by transforming the pressure records using linear 
theory. 

A mean incident wave angle (a) and a representative fre- 


quency (f,) are defined for input into the models such that 


49 





using a and "5 Poscatleurace Ene ytotal radiation stress Sey? 


the same value is obtained as that measured by the pressure 
sensor array. 
The radiation stress ($7) in equation (4.4) can be ex- 


pressed applying linear wave theory for long-crested waves 


oxy t? = EB(fi@n(f) sin a cos a (46) 


The total radiation stress (S07) is calculated using 


S = f Syy (Ff) dé (4.9) 


if 
y 


in (4.8) for monochromatic (f = 5) waves approaching ata 


Then, the calculated Sey is equated to its equivalent S.. 


mean angle a 


Ay ale _ — 
> big e 
E AM sin a(t) cos «a ( =! (A. 004 


where aye is the total energy. Applying linear theory rela- 


tionships a mean incident wave angle a is defined as: 


a(£.) = + sin?+ js, (Ae) 
P a [eT (ene) | 


where the representative frequency ee is the peak frequency 
of the energy spectrum. Since the radiation stress is calcu- 


lated using the pressure gradients, the kinetic energy calculated 
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using the pressure gradient is introduced for consistency 
of the analysis. It is assumed the kinetic energy Ey 1s 


half the total energy. Then equation (4.11) yields 


gl 
nle 92 >sin~ are (in 


p B,ntt) 
The kinetic energy is calculated using the pressure gradient 


Semponent spectrum 


f) = fH(£) [1G (£) + Gy (£)) (45435) 


aX oY 


where lH (£) |? is the linear wave transformation function 
given by 
2 ae sinh (2kh) 
lH(#)|° = 22 —ss (4.14) 


Aku“ cosh*k(h +z_) 


The mean incident wave angle (a) at the pressure sensor's 
location was then refracted shoreward to the 5 m contour (at 
the shallow water limit). A refraction calculation was per- 
formed over a bathymetric chart of the experimental site that 


shows the beach has straight and parallel contours shoreward 


of 5 m depth (Fig. 4.4). Applying linear refraction to the 
measured mean angle (Ae eat the pressure sensor location 
and peak frequency (EE = 0.07 Hz), waves reach the 5 m contour 
with a refracted angle O. = 9.0° where the calculated refrac- 


tion coefficient K equals 0.985. 
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. SANTA SAREFFPA 


DIMENSIONS IN METERS 
REFRACTION __ 2'A GRAM 


Figure (4.4). Refraction diagram for Santa Barbara 
(17 February, 1980). 
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The ae is calculated using linear refraction theory: 








C 
Ht ms J gm ee oon ye 
rere ncoE ae 
rms (8m) g 
“Sp ne 
a haa Time (KL) (4215) 


sl 


Giving a refracted H equivalent to .85 m at water depth 


rms 


De = 5 m, which is used as the initial wave height for the 


model. Refracted wave rays were extended to get the angle 


of approach at the breaker line calculated to give Cy = 602. 


Mies breaker line was defined at maximum eee Ch was used 


to calculate the longshore sand transport using the wave 
power expression. 
Table II gives various calculated parameters of Sey! a, 


uv 
Ef — at different times before, after and during the 
17th and 18th of February. Orientation of the coordinates 
deads to a negative sign for both Sh and a, indicating the 


waves are coming from the west. 
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TABLE. if 


Calculated Parameters at West Pressure Sensor 
(During 17-18 February) 


Date | time |depth tide S = H 
XY 5 p rms 
(m) (cm) (Eig enn) (Hz) (cm) 
5 
a7 dl -5.8 65.6 
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V. LONGSHORE SEDIMENT TRANSPORT MODEL 


A littoral sand transport model is developed based on 
the energetics approach expressed by Thornton (1973) in which 
the energy transformation is described by the random wave 
model of Thornton and Guza (1983a) coupled with the longshore 
Seerent model of Thornton and Guza (1983b). The input to 
the predictive model are the offshore rms wave height, mean 
direction and mean period. Both an analytical model for plane 
sloping beaches and a numerical model for arbitrarily bottom 


profile are presented. 


A. DEVELOPMENT OF SAND TRANSPORT MODEL 
Assuming stationary wave conditions with a beach of 


straight and parallel contours, the gradient of the wave 
dECg 
dx 


ensemble averaged dissipation due to both wave breaking 





power in the on-offshore direction is balanced by the 


<€,> and bottom friction <e¢> (Thornton and Guza, 1983) such 


that 


dECg, 
dx b ig 


The wave power is expressed using the linear wave theory 


relationships: 
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a, Ee 2a _ 
Cg, = 7 [l — Raa COs tS ta.) 


a 


where H__. is the root mean square wave height, a is the 
mean wave direction and k is the wavenumber associated with 


peak frequency Eo: 

Dissipation due to wave breaking <E,> can be modeled 
after a periodic bore, applying conservation of mass and momen- 
tum at regions of uniform flow up-and-downstream of the bore 
as shown in Fig. (5.1). The average rate of energy dissipa- 


tion per unit area for single breaking wave is calculated 


(Stoker, 1957): 


(hea hs) 4 
L 2 i ieee H 
ar ©. Sl) = a eet Ol oe ie ©? (Oey (5.4) 


€ 
bore 


where the wave height H is measured as the maximum to mini- 
mum of the bore, Q is the volume discharge of flow per unit 
area across the bore described by simple linear bore theory 
after Hwang and Divoky (1972), and B is a breaker coefficient 
of O(1). The coefficient B is an unspecified parameter in 
the model, which is determined from the data. 

The Rayleigh distribution is assumed to describe the 
waves everywhere with the implied assumptions that the waves 
are very narrow band such that they can be described by a 
single frequency f and mean direction a. Thornton and Guza 
(1983) predict the fraction waves which are breaking at 


each location. Using this information, the average rate of energy 
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Figure (5.1). Periodic bore used to describe 
spilling breakers. 


a 


CS aaa 





dissipation due to breaking Soe 1s given by 


3 
- Bet 
4605 ¢= 3v7 H/ F C52 5) 


=a 2 So 
b 16 ae rms 


where y 1S a coefficient relating eee to h (taken to be 0.44 


by Thornton and Guza, 1983). The average frictional dissi- 


pation <e€.> at the bottom boundary layer for the ensemble 


f£ 
(all waves) calculated using Rayleigh probability distribution 





is: 
le) Ce 21 OE ; aye 3 
oa (Srna ae ee! 
l6vn 
where C. is the friction coefficient. Thornton and Guza 


1g 


(1983a) point out that <e.> is relatively more important for 


fe 
gently sloping beaches and for low frequency waves, i.e., the 
shallower the depth the greater the frictional dissipation. 
They neglect the frictional dissipation as it has a small con- 
tribution compared with <€,>, as will be done here also. 


Po seleurtiicg equations (5.2), (5.3), (5.5), and (5.6) 


Paco equation (5.1), yields: 


a. 
3¥y 8 GF, OB 7 


ainsi Z. —_ | 
dx 8° 9 Brms Cy a 7 LG a Hyms '?-!? 
ae can be evaluated analytically, assuming small inci- 
dent wave angle (cos a = 1) and limiting the analysis to 


ss =), which gives 


ty 5° 


Shallow water ( 





Z leo 7 0 Bec) 4 1 2 -1/5 
a = a h [kee la 3374 Eye. ’ CSc.) 
O Yo 
where 
eee lye 
Can Ho hs’ 
and 
— 2g (gy) 142 a tan 8B 
5 IRS. a 3 


Bef 
P 


The value of tan 8 is the beach slope and the subscript "o" 
refers to the input conditions at the outer shallow water 
depth limit. When the water depth gets very shallow, (5.8) 
becomes 
_ Base 7 ash 0 (5.9) 
which indicates that the wave height in the inner surf zone 
is related to the depth and independent of the initial condi- 
tions in deeper water. 

With the assumptions of straight and parallel depth contours 
and stationary wave conditions, the longshore current dis- 
tribution across the surf zone can be estimated using the 


longshore momentum flux balance (Thornton and Guza, 1983) 


dS 
X 


ar ys (57, 10) 


Sy, 





where T, is the bottom shear stress and Sey mow the Lracgclation 


b 
stress defined by Longuet-Higgins and Stewart (1962) (for 


linear theory) as 


S =e C COS. a east 
XY g . 


From linear wave refraction 


Sing o2 oe 
—— ao OO constant 
O 
Thus, 
dS, sin a, fg 
—— = or lac i rs Cos a | Ga 


The term in brackets in the R.H.S. is defined in (5.1) 





as <E)>- Therefore (5.11) reduces to: 
aS sin oc 
Vaew 2° 3 
Ee = c- Sea = Th (Se) 


The average alongshore bottom shear stress Tt, is defined by 


b 


Longuet-Higgins (1970), assuming the quadratic law 
+ 
(V + Vie) (oes) 


where a, 1s the total instantaneous velocity due to waves 
and current, V is the longshore current, ve is the longshore 


component of the oscillating wave velocity and Ce is the bed 
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shear stress coefficient. The bar over equation (5.13) indi- 
cates averaging in time over one wave period. Assuming a 
small angle of wave incidence (a < 15°), and a weak current 
case following Liu Ph. L-F and R. Dalrymple (1978), equation 


(5.13) simplifies to 


[ S 6 Ge hy (S24) 


where u is the mean cross-shore wave velocity. 
Ju. is specified using the shallow water wave theory 
relationships and Rayleigh wave height distribution (Thornton 


and Guza, 1983) such that: 


= glee AZ 
= aay! H ems (5-15) 


Biles cimeic Uiqmm@erenl 1) pemeSa 12) 58e(5.14)40(5.15) intoew(5.10) 


and solving for the mean longshore current velocity, yields: 





Sin a 
y= A = = O 46. 
p cy a O 
3 B> fc gi’? Sin 4° 
a ee aa ren} 
4 .. ea Cy n2/2 


The longshore sand transport distribution is obtained 


by substituting the expressions for the change in energy flux 


JdEC 
hy equation (5.1) and the longshore current equation 


feo) into equation (2.11) to give 
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B sin a WAP Bye 
Y ee ons <e,> 6517) 
g(1-5-) 14, f O 


Pelosi reuitrng (9.7), (5.15) and (5.16) in equation (5.17) 


yields: 
“a B l sin aS fe B? Ho lug 
ee (ete) 
g {1 ===) 7 <I9) O ia 


Ss 


Equation (5.18) has the observed characteristics that outside 
meme surt zone as h reaches hoe q. (x) gets small and when h 
goes to zero at the shoreline, q. (x) approaches zero. A 
maximum transport occurs between these two boundaries. 

For a planar sloping beach, the analytical solution for 
ec (5.8), which is limited to shallow water, is substituted 
maton ( 5.18) to obtain the cross-shore transport distribution. 
Heieearoitrary bottom profiles (but straight and parallel depth 
contours), a simple forward differencing scheme is used 
to numerically integrate equation (5.1) for Has? The 
calculated See Walnvewcesubstitrubed Into (5.18). The energy 
flux balance is also numerically integrated from initial 


offshore distance (at depth ho) to the shoreline with a 


Spatial interval Ax 


ECg,li4, = £ oe |= ax{<e,>|, - cee lil (Eyes) 


a eee 2 
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Equation (5.19) starts at the initial offshore distance 


; fi Ot are given. The model calculates 


the forward step by subtracting the total dissipation from 


(1 = 1) where Hos A, 
the calculated energy flux evaluated from the previous step. 
The model predicts eS each step, uSing linear theory rela- 
tionships (5.2) and (5.3) until i approaches n at the shore- 
line. Then, the predicted He values are used to calculate 
the longshore current using (5.16). The sand transport dis- 
tribution is evaluated using (5.18). 

Integrating qe in (5.18) over limits of offshore distances 


for the study zone, the volume transport rate (Q ) 1s calcu- 


lated as 
Q = jf qd, dx (Seen) 
where XO is the offshore distance at which the initial condi- 


tions of waveheight (Ho) and water depth (ho) are given. 


The coefficient Be in (5.18) is used to calibrate the model. 
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VI. MODEL RESULTS AND COMPARISON WITH PROTOTYPE DATA 


Pre MODEL RESULTS FOR A AND V 
rms 


Field data for the Santa Barbara experiment are utilized 


to calculate the parameters S_,a_, f. and H which in 
ye EO p rms 


turn are used as input for both the expressions (5.8) and 


(5.19). Model coefficients y, B and Cc. were determined from 


H 
the data. The ratio, y = = was found to be 0.44 from the 





slope of the FA vs. sae curve, Fig. (6.1), in the inner surf 
zone. This agrees with the value reported earlier by Thornton 
and Guza (1983) for the Torrey Pines Beach experiment. 

ens values inferred from the pressure and velocity 
measurements at time 1800 on the 17th of February are compared 
with the wave transformation model for various breaker coeffi- 
cients, B. Using an iterative scheme, based on the least 
squares error criterion, B equal to 1.1 gave the best fit to 
a data (Fig. 6.2a). To check on the generality of the 
model, the waveheight transformation model was run for 5 Feb- 
ruary using B= 1.1. It also gave the best model fit to 
the data (Fig. 6.2b). 


An optimal bed shear stress coefficient, c was calculated 


Bed 
comparing the longshore current model with the measured long- 
shore current values, V. Since there are only three longshore 
current measurements available on 17 February due to the des- 


Eruction of instruments by the storm, the data from 5 Feb- 


ruary were again used. With B= 1.1, Ce was varied to 
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y = 0.44 in the inner surf zone. 
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Figure (6.2). Model H as a function of offshore distance 


compared with measurements indicated by x. The 
model is applied to data; (a) on 17 February; 
(b) on 5 February, to obtain B which was 
feOMnGwtO fea.) yin both days. 
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obtain the best fit to the measured V. A value of Ce = 6 8S 


gave the best fit (Fig. 6.3b). 


pee MODEL RESULTS FOR CROSS=SHORE SAND TRANSPORT DISTRIBUTION 
The empirical coefficient Be in equation (5.18) is used to 

calibrate the model and will be explained later. In all the 

calculations and figures, the offshore distance (x) is 

measured starting from the experimental baseline at x = 0. 

The average net transport between profiles on 17 and 18 

February are to be determined. A weighted time difference 

of 25 hours was found between the two mean profile surveys. 


The model predictions of H V and ae at each x-location 


Ems 
depend mainly on the water depth (h). Therefore, two cases 
are considered to test the sensitivity of the models to the 
variations in the water depth. In the first case, the water 
depth does not vary with time. In the second case, the water 
depth varies by hourly time increments according to tidal 
Variation. 
1. Model Results for Constant Water Depth 

In this case study, the water depths are held constant 

aeecelative MSL. The numerical model utilizes the actual 


bottom and predicts H V and qe: The model results are 


rms’ 
summarized in Fig. 6.4 indicating the capability of the 
model to predict the increase in waveheight due to shoaling 
and the decrease due to breaking processes. The longshore 


sand transport and the longshore current have a similar dis- 


tribution, which implies that the sand transport may primarily 


on 


80 


V (cm/sec) 


80 


60 


40 


V (cm/sec) 


meaure (6.3). 


50 x 100 
X (m) 


150 200 


> 


50 100 150 = * «200 


3 X (m) 


Longshore current (68-minute average) as a 
function of offshore distance compared with 
measurements indicated by x. The model is 
applied to data on; (a) 17 Feb., showing only 

a few instruments working; (b) 5 Feb., for which 
Ce = 0.005 was found to give the best fit. 


68 





be governed by the longshore current (Fig. 6.4). The influ- 
ence of rapid depth changes is reflected in the small pertur- 
bation of the V and qe curves at x = 65 m, which is associated 
with a local flatness of the beach slope, as can be seen in 
the "mean beach profile," Fig. 4.1. The maximum waveheight 
is used to define a "mean breaker line" location, Xpe jRepe 
discussion purposes and normalization of the results. The 
analytical model is applied to a plane sloping beach. A slope 
of 0.037 was found as the best fit to the "mean beach" pro- 
file using Fig. 4.1. The normalized distributions resulting 
from the analytical model are summarized in Fig. 6.5, giving 
the maximum sand transport at x/X, = OL G0 se whtel is Simila: 
to the results of Komar (1977a). 

Comparison between the sand transport distributions using 
the numerical and analytical models is given in Fig. 6.6. 

The sand transport peak of the numerical model, which 
uses the actual bathymetry, is shifted offshore compared 
with the analytical distribution using a planar sloping beach. 
It may be seen in equation (5.18) that qe is proportional 
to =), which means qe increases as h decreases for the 
same location. The actual bottom profile (Fig. (4.1)) used 
in the numerical model is deeper near the breaker line com- 
pared with the plane sloping beach used in the analytical 
model; therefore, the calculated maximum q. in the numerical 


model is shifted offshore relative to the analytical model. 
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The models predicting H__ 5, V and qe were run for 
the 25 hours using a one hour time step and correcting the 
depths for the tidal elevations at each time step. The 
resulting 25 sand transport distributions were averaged over 
the 25 hours to give a predicted mean longshore sand trans- 
port distribution between profile surveys on 17 and 18 
February. 

Mean sand transport distributions including tidal 
variation of both the analytical and numerical models are 
given in Fig. 6.7 which shows a good agreement with a peak 
shift towards offshore for the numerical model following the 
previous discussion. When considering the tides, the averaging 
causes the distribution to spread onshore to x = 30 m with 
less ce than the first case for the same coefficients. 

In both cases the maximum longshore sand transport occurs at 
almost the same location as the longshore current which indi- 
cates again that the sand transport is primarily driven by 
the longshore component in both models. The sediment trans- 
port distribution predicted by the numerical model is not 
smooth, unlike the analytical curve. This is because of the 


use of the actual depth profile. 


fee ODEL RESULTS AND COMPARISON WITH PROTOTYPE FIELD DATA 
The net longshore transport rate distribution for Santa 


Barbara Beach is calculated using the profile changes during 


TS 
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Figure (6.7). Average sand transport distribution over 
the daily tidal cycle calculated using 
variable depth; for analytical model 
( ) and numerical model (----). 
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17 and 18 February. Most of the estimated sediment trans- 
port occurs at a distance 20 m from the baseline indicating 
a Significant beach face erosion, which is beyond the scope 
mecie present study. 

Since the volume sand transport rate calculated from the 
data is the net volume change rate in the beach profile, and 
does not represent the total longshore sand transport as the 
model predicts, care must be taken in the interpretation of 
the comparisons. 

A dimensionless sand transport distribution for the field 
data 1s compared with results from the numerical model in 
Fig. (6.8). The measured maximum net longshore sand trans- 
port occurs on the beach face at x = 20 m. The hypothesis that 
the predicted longshore transport distribution is proportional to 
the net sediment transport distribution, calculated from 
the profile changes does not appear to be the case. Signifi- 
cant profile changes occur at the beach face which are not 
included in the model calculations. It is also possible 
that the sand eroded from the beach face was moved in the 
offshore direction outside 100 m and was not included in the 
measured profiles, which would also invalidate the original 


hypothesis. 


fee init EMPIRICAL COEFFICIENT Be AND LONGSHORE SAND TRANSPORT 
The empirical coefficient Bor as described by Thornton 


(1973), has a magnitude of order (1) and is equal to 


(es 
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Figure (6.8). Normalized sand transport distribution 
inferred from profile changes vs. 
numerical model results with variable 
depth. 
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i (6.1) 





where e is the transport efficiency of the available wave 
power and fy is the dynamic friction factor associated with 


sediment transport. The factor fy 1s equal to: 


- sediment shear stress 
h normal stress at the bottom ’ 
while b is a proportionality coefficient between the ratio of 
sand velocities (=) and the sguare root of the water particle 
maLocity ratio a. 

VW 

Because the data do not provide information about each 
individual parameter, all parameters are combined 
Te" Bo, which is determined empirically by model calibration. 

The total longshore sediment transport rate (Q,) ner 


(m>/day) is calculated by integrating the area under the sand 


meansport distribution 


The calculation based on profile changes does not give the 
total longshore sand transport but gives the net transport 
which is much less in this case than the longshore transport 
predicted by the model. Therefore, equation (2.5) suggested 
by CERC (1977) is used as a reference to calibrate the long- 


shore sand transport model. 


aa 





The calculations of the longshore sand transport rate Q, 
require the wave condition at the breaker line ee AL 
and Oe Field data give es equal to 0.98 m, Fig. (6.2a) 
and Ay to be 2.96 m knowing the beach slope besides the angle 
Oy, equals six degrees from the refraction diagram, Fig. (4.4). 
Q, calculated using equation (2.6) gives 4600 m°/day, which 
requires B to be 0.18. 

A summary of the results if given in Table III. The 
difference between the two cases of depth (constant and 
varying) for the analytical model, was found to be insignifi- 
cant. A similar result was found for the numerical model. 
The calculated Qy using the numerical model is less than Q, 
calculated from the analytical model, because the numerical 
model used the actual bottom profile. The numerical and 
the analytical models gave essentially the same maximum 
cross-shore transport distribution Cees for the same case 
of either constant or varying depth, Fig. (6.6) and Fig. 
or?) . 

Johnson (1953) estimated the longshore sand transport 
at Santa Barbara Beach to be 220,000 m°/year, Or vabOue O00 
m /day on the average. The Qo predicted by the model gives 
almost eight times the above estimated average daily rate. 


The high value of Qo predicted by the model reflects the 


storm conditions, which caused significant longshore transport. 
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VII. SUMMARY AND CONCLUSIONS 


The principle objective of this study is to predict 
the longshore sediment transport distribution. Special 
emphasis was placed on testing an analytical and numerical 
model against field data acquired from an experiment at 
Leadbetter Beach, Santa Barbara during the intensive storm 
period 17 to 18 February, 1980. 

Field data show a strong longshore current coupled with 
Significant changes in the beach profile with no evidence 
of on-offshore sand movement (bars or troughs) along the 
profiles. Therefore, it is hypothesized that the changes 
in the profiles are due to the variations of sand transported 
alongshore by the observed longshore current, and that the 
on-offshore sand transport is negligible. The net sand 
transport calculated from the profile changes is expected 
to have a distribution similar to the cross-shore sand 
EitrisoOrt distribution. 


Initial values for the models of He + Ot and f are 
O 


obtained from the field data. Model coefficients y = 0.44, 
Bee -L.l and Ce. = 0.005 are determined by model fitting. The 
developed model is capable of predicting the increase in wave 
height due to shoaling and the decrease after wave breaking. 


The model predicts the longshore current (V) and the longshore 


sand transport (q.) at each location. A maximum q. is 
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predicted at distance x = 0.8 x, which agrees with Komar 
(1977a). The experimental transport coefficient, Bes ist Que 
mompe O28 anstead of 0-08 by Thornton (1973). The maximum 
longshore current, Vee, occurs at approximately the same 


location as (q_) indicating that the sand transport 


s‘’ max’ 
distribution is primarily determined by the longshore current 
fers tribution. 

The net change in the beach profiles used to infer the 
net sand transport distribution, shows qe maximum at the 
beach face. Since the sand transport model does not consider 
the swash region, the net profile change distribution does 
not agree well with the predicted cross~-shore distribution, 
and the original hypothesis does not appear to be valid. As 
a result, sand transport mechanics in the swash zone appear 
to be important in determining why the maximum sand trans- 
port occurs at the beach face. The swash is driven primarily 
by infragravity waves (long waves) due to surf beat and edge 
waves which have not been included in the surf zone models 
where only the depth-limited sea-swell frequency band of 
waves are considered. Therefore, further studies need to 
model the swash zone in order to describe the onshore water 


wave boundaries in amore realistic manner. 
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